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This report is documented in two volumes. Volume I describes modifica¬ 
tions made to the Mark IV Supersonic-Hypersonic Arbitrary-Body computer pro¬ 
gram, and Volume II documents two computer graphics codes y£ed/to validate 
Mark IV geometries. 
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SECTION 1 


INTRODUCTION 


This report describes new developments related to the Mark IV 
Supersonic-Hypersonic Arbitrary-Body Program (Reference 1), a FORTRAN 
computer code employing design methods for computing the aerodynamic 
characteristics of complex configurations. The effort is documented in 
two volumes. Volume I describes improvements made to the Mark IV 
program itself, and Volume II documents two computer graphics codes 
used to validate Mark IV geometries. Contained in Volume II (Reference 2 ) 
is a detailed description of TEKPIC, a new interactive graphics code 
that allows the user to examine many orientations of a given vehicle in 
a minimal amount of time. However, the TEKPIC program uses only an 
approximate method for displaying the visible lines of the geometry. A 
second graphics code, HIDDEN, requiring significantly more central pro¬ 
cessor (CP) time and core memory than TEKPIC, but capable of removing 
hidden lines, is used to complement the TEKPIC program. 

The primary objective of the work documented in Volume I was to 
incorporate into the Mark IV program a method for tracing inviscid 
surface streamlines over arbitrary geometries. The Newtonian streamline 
method, which provides the flow direction at any point on the surface 
given only the freestream velocity (vector) and the local surface 
outward normal, is employed in the modified Mark IV program. 

Although the basis of the streamline method is simple, many diffi¬ 
cult problems arise when attempting to generalize the approach to 
arbitrary geometries. 

Normally, the inviscid flow properties (including the direction 
cosines of the surface velocity) are known only at specific points on 
the surface. Since streamlines do not generally pass through these 
points, a method must be available for interpolating between the points. 

The interpolation procedure becomes quite complex for arbitrary configurations. 
Another problem encountered is that of locating the origins for the 


1 




streamlines. Streamlines may emanate from blunted stagnation regions in 
which the surface outward normal is aligned with freestream velocity, or 
they may originate from the leading edges of wings, canards, nacelles, etc. 
in which the surface outward normal is not aligned with the freestream 
velocity. Given only the orientation of the outward normal with respect 
to the freestream vector, a consistent method can be developed for 
locating stagnation region origins but not the latter origins. 

An initial attempt to develop a general-body streamline method was 
made by the authors of the Mark IV program. A surface spline technique 
(Reference 3) was applied to certain user-specified regions of the geometry 
as a means of interpolating for the surface flow properties. (The term 
region, as used in the original Mark IV program, is defined as a collec¬ 
tion of geometry Panels; see Reference 1, Volume I for a complete defini¬ 
tion of Panels, Sections, and Elements.) The surface fitting method has 
the disadvantage that it is not a parametric spline. That is, the surface 
fit of a given flow quantity within a region requires an appropriate choice 
of two of the three independent variables to ensure that the dependent 
variable is not multivalued. Each region must therefore be fit indepen¬ 
dently of other regions, and positional continuity between spline regions 
is not guaranteed. This poses a problem when attempting to trace stream¬ 
lines from region to region. The latest release of the Mark IV program 
contained no means of tracing streamlines across region boundaries. There¬ 
fore, running lengths along the streamlines could not be accurately pre¬ 
dicted, and no method could be implemented for locating streamline origins. 

Although the surface spline used in the Mark IV program has several 
disadvantages, the incorporation of other more sophisticated spline 
methods would be beyond the scope of the present effort. Furthermore, the 
use of other spline methods would require significant modification of the 
geometry input procedure which was not to be altered. However, in the new 
streamline method, improvements were made to the manner in which the 
nonparametric spline method was used. Rather than surface fitting 
quantities over large regions of the geometry, which may lead to erroneous 
interpolations, the surface spline is applied only to each geometry Panel. 
Not only does this approach improve the accuracy of the interpolations, 
but it greatly simplifies the tracing of streamlines from one spline 
region (Panel) to another. (The term region, as used in the remainder of 
this report is synonymous with the surface of one Panel.) 


2 


The new streamline method also has the ability to locate the origins 
of the streamlines. This is accomplished by distributing starting points 
along the aftmost boundaries of the geometry, and tracing the streamlines 
forward, against the flow direction, until the appropriate origins are 
located. After all streamlines are computed, the distribution is 
examined to ensure that it is sufficiently distributed over the geometry 
for subsequent boundary layer calculations. If some surface areas are 
void of streamlines, new starting points are strategically positioned and 
more streamlines are traced. 

The user may wish to know surface property information (including 
viscous related parameters) only at specific points on the body. An 
option has therefore been included which allows starting points to be 
input by the user. If starting points are input, the user must request 
whether the integration is to be performed with or against the direction 
of flow. If the streamlines are integrated in the direction of flow, it 
is assumed that the starting point is also an origin. 

The primary purpose of calculating inviscid surface streamlines is 
to provide paths along which boundary layer methods may be applied. A 
secondary objective of the effort described in this report was to ensure 
that the integral boundary layer methods employed in the original version 
of the Mark IV program were compatible with the new streamline method. 

The integral methods (References 4 and 5) were originally coded in a 
separate computer program by McNally (Reference 6) who was mainly concerned 
with boundary layers in shock-free flowfields. In the coding of the 
equations, isentropic conditions were assumed to exist along all streamlines 
extending from the freestream to the geometry surface. Therefore, the 
McNally coding required modification for use in the Mark IV program. 

However, several discrepancies were discovered in the Mark IV coding of 
the integral methods. Many of the boundary layer edge quantities appearing 
in the integral equations were still based on freestream conditions instead 
of local conditions, e.g. p Q , p . Therefore, an additional objective of 
the effort described in this report was to correct the coding of the 
integral methods. 
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One constraint placed upon all modifications to the Mark IV program 
was that the overall operation of the program was not to be affected. 

Only those parts of the program directly involved with either the 
streamline tracing or the integral boundary layer methods were changed. 

No modifications were made to the geometry package, the inviscid aero¬ 
dynamics methods, the shielding analysis, or the special routines. As 
a consequence, unless streamline or viscous calculations are desired, 
the user follows the input data formats exactly as they are described in 
Volune I of the original Mark IV docimentation. 

The streamline method and the viscous analyses are accessed from 
the AERO executive routine in the same manner as described in the original 
documentation. Section II of the present report describes the input data 
required to use the new streamline method and the modified viscous 
methods option. The information in this section is intended to completely 
replace the input data instructions for both the Surface Streamline 
Option (pp. 92-98 in Volume I of Reference 1) and the Viscous Program 

Option (pp. 114-128). However, the format of the Mark III Program 
(Reference 7) Skin Friction Element Data Cards used in the Mark IV 
program, has not been altered. The only other data formats changed 
in the Mark IV program were those associated with the Surface Data 
Transfer Option, used only when streamline calculations are desired. 

A detailed description of the theory and the algorithms used in the 
new streamline method is given in Section III of this report. A 
relatively brief description of the integral boundary layer equations used 
in the Mark IV code is also presented in Section III. The theory and the 
coding of the integral methods are well-docimented in McNally’s report 
(Reference 6). Section IV of this report contains general information 
for those wishing to modify either the streamline method or the viscous 
methods. 

The modified Mark IV program is operational on the CDC CYBER 750 
computer maintained by the ASD Computer Center, Wright-Patterson AFB, 

Ohio. 


4 



SECTION II 


USER'S GUIDE 


The user's guide presented here is intended to serve only as a 
supplement to the original user's manual of the Mark IV program 
(Reference 1, Volume I). Only those changes that affect either the 
input data formats or the operation of the original version of the 
code are documented in the present report. Unless the user is 
interested in tracing streamlines or in computing viscous effects, the 
user may rely solely on the original user's manual. 

The general organization of the modified Mark IV program, shown 
in Figure 1, is identical to that of the original version of the code. 
In the normal operation of the program the geometry package is called 
first to generate and save on local storage unit 4 the quadrilateral 
element data. The geometry data on unit 4 (a random access or mass 
storage unit) includes the four corner points, the three components of 
the outward normals, the areas, and the centroids of all elements 
comprising the vehicle geometry. 
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Once the geometry data have been saved on unit 4 the user may proceed 
with the aerodynamic calculations. The aerodynamics program, AERO, con¬ 
sists of six separate components, and each component is responsible for 
performing a certain analysis. The AERO program requires the user to 
specify a sequence of integers identifying the particular analyses to be 
employed and the order in which they are to be used. Technically, the 
analyses may be called in any sequence, and each analysis may be called 
any number of times. In practice, however, some analyses depend on re¬ 
sults generated by other options in AERO. For example, the Special Rou¬ 
tines option, responsible for summing the forces and moments of the 
various vehicle components, requires that the Inviscid Pressures analysis 
be performed first. 

Prior to accessing the Streamline option it is necessary that a 
specific sequence of calls be made to other AERO options. The calculation 
of surface streamlines requires the prior computation of surface flow data. 
Two methods may be used to generate the surface data. The data may be 
hand-loaded using the Surface Data Transfer sub-option of the Flow Field 
Analysis option, or the data may be generated using the Inviscid Pressures 
option. The hand-loaded data are stored on unit 10, a random access unit 
containing all flowfield information generated by the Flow Field Analysis 
option. The formats of the surface data on unit 10 are compatible with 
the Streamline option, but the data generated by the Inviscid Pressures 
option are stored on unit 4 in formats not compatible with the Streamline 
option. Therefore, if the surface data are generated by the Inviscid 
Pressure option, the data must be transferred from unit 4 to unit 10 be¬ 
fore the Streamline option may be called. The data transfer is accom¬ 
plished by the Surface Data Transfer sub-option of the Flow Field Analysis 
option, the same option used to hand-load surface data. The surface data 
on unit 10 may consist of a combination of data generated by the Inviscid 
Pressures option and hand-loaded data. Once all the surface data are 
stored on unit 10, the Streamline option may be accessed. 

The utilization of the Viscous Methods option also requires that 
other analyses be performed first. The Viscous Methods option contains 
two approaches for computing boundary layer effects. The first method 
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is the Mark III Element Skin Friction method which requires that the 
Inviscid Pressures option be called previously. Although the Mark III 
Element Skin Friction analysis is fairly crude, it is the only method in 
the Mark IV program capable of estimating the contribution of skin 
friction to the total vehicle forces and moments. The second approach 
available in the Viscous Methods option applies boundary layer methods 
along the inviscid streamlines generated previously in the Streamline 
option. Although the approach is more sophisticated than the Mark III 
Element Skin Friction method, it is not capable of predicting the contri¬ 
bution of skin friction to the total vehicle forces and moments. The 
latter method should be used only when detailed boundary layer informa¬ 
tion is desired (momentum thickness, displacement thickness, etc.), or 
when heat transfer prediction 0 at specific points on the body are 
required. 

Shown in Figure 2 are the three possible user-specified sequences 
of AERO calls to the various options that may be used to generate stream¬ 
lines and to subsequently make integral boundary layer calculations along 
the streamlines. The particular sequence of calls used is specified on 
the Aero Flag Card described on p. 69 in Reference 1, Volume I and are 
read by the AERO program. Following the Aero Flag Card are the Flight 
Condition Card, the Reference Dimension Card, and the ot-6 Cards which are 
also read by the AERO program. The remainder of the data to be prepared 
depend upon the specific AERO options selected. For example, if the first 
option is the Inviscid Pressure option, the user would turn to the Pres¬ 
sure Calculation Program Input Data on p. 101 in Reference 1, Volume I. 
After all data for the Pressure Calculation Program are prepared, the 
data formats for the second option are located in the user's manual, and 
so on. However, if the option selected is the Streamline Analysis option, 
the Viscous Methods option, or the Flow Field option in which surface data 
are to be hand-loaded or transferred from unit 4 to unit 10, the user must 
refer to the following descriptions of the data preparation. For optimal 
use of the Streamline and Viscous options it is suggested that the user 
!■■ coo me familiar with the theory of the methods presented in Section III. 
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1. SURFACE DATA TRANSFER OPTION 


The data formats given for the Surface Data Transfer option in Ref¬ 
erence 1, Volume I (pp. 89-91) may not be used with the modified Mark IV 
program. Instead, the user should refer to the data formats presented here. 

This option is one of several in the Flow Field Analysis program, and 
is used to transfer surface data to unit 10 from either the input unit 
(hand-loaded data) or unit 4 (data generated previously by the Inviscid 
Pressures option). The data are placed on unit 10 in formats acceptable 
by the streamline option. The Surface Data Transfer option is exercised 
only if IDTYP(l) = 2 on the Region Directory Table Card, Reference 1, 

Volume I (pp. 76) read by program FLOW, the executive program for the Flow 
Field Analysis option. This option may also be used to read and print out 
surface data previously stored on unit 10 if IRW = 1 on the Region Directory 
Table Card. 

As described in Reference 1, all flow field data generated or read 
by the Flow Field Analysis option are stored on unit 10, the flow field 
unit. Each of the analyses in the Flow Field Analysis option, including 
the surface data transfer, may be applied to any user-specified collection 
of geometry Panels. For a given angle of attack and set of freestream 
conditions, the results of each flow field analysis are stored in a 
"flow region." The region number, IREG, associated with a particular 
analysis is identified by the user on the Region Directory Table Card. 
Surface data transferred from unit 10 must be stored in flow region 1 
(IREG =1), and hand-loaded surface data must be placed in flow region 2 
(IREG = 2). If streamline calculations are desired, the surface data in 
regions 1 and 2, taken collectively, must uniquely define the surface 
flow quantities over the surface of the entire vehicle. 

When transferring surface data to unit 10, the information is further 
divided into "subregions." A subregion is simply defined as the collection 
of surface flow properties over one geometry Panel. In the original ver¬ 
sion of the Mark IV program, the surface properties associated with sev¬ 
eral Panels could be grouped together into one subregion. However, since 
the surface properties associated with each subregion are fit with a sur¬ 
face spline (for interpolation purposes) in the Streamline option, erroneous 
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interpolations may result if a subregion contains the properties of more 
than one Panel. Therefore, the number of subregions to be loaded, NSREG, 
as specified by the user on the Flow Field Control Card below, simply 
corresponds to the number of Panels whose surface property data will be 
placed on unit 10. Only one Flow Field Control Card is read for each call 
to the Surface Data Transfer option. 


Flow Field Control Card (312) 

This card is input only if IDTYP (1) = 2 and IRW = 0. 


Column 

Code 

Routine 

Format 

Explanation 

1-2 

NSREG 

FFSURF 

12 

Total number of subregions (geometry panels) to 
be loaded from unit 4 to unit 10 or hand loaded 
(assumed at least = 1). 

3-4 

KTRNSF 

FFSURF 

12 

Type of surface data transfer. 

= 0 Data will be read from unit 4 and placed 
on unit 10. (Next card read will be the 
Surface Data Panel Selection card). 

t 0 Data will be hand loaded and stored on 
unit 10. (Next cards read will be the 

Hand Loaded Surface Property Data cards). 

5-6 

IPRINT 

FFSURF 

12 

Print flag. 

= 0 Do not print surface data. 

= 1 Print. 


Hand-Loaded Surface Property Data Cards 

These cards are used to hand load surface property data directly to the flow 
field unit 10 and are used only if KTRNSF f 0, otherwise skip this section 
and go to Surface Data Panel Selection Cards, pp. 12. The first card identi¬ 
fies the gometry Panel with which the data are to be associated and specifies 
the number of data points (sets of Surface Data Coordinate and Surface Data 
Property Cards) to be read. The number of sets of Surface Data Panel 
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Hand-Loaded Surface Property Data Cards (Continued) 

Identification, Surface Data Coordinate, and Surface Data Property Cards 
must equal NSREG on the Flow Field Control Card. 


Surface Data Panel Identification Cards (215) 


Column 

Code 

Routine 

Format 

Explanation 

1-5 

IPANL(I) 

FFSURF 

15 

Panel number on unit 4 with which the data are 
to be associated. 

6-10 

NPTS 

FFSURF 

15 

The number of data points to be read for this 
subregion. 

Surface 

Data Coordinate Card (6F10.0) 

1-10 

DATA(l) 

FFSURF 

F10.0 

X-coordinate of the surface data point. 

11-20 

DATA(2) 

FFSURF 

F10.0 

Y-coordinate of the surface data point. 

21-30 

DATA(3) 

FFSURF 

F10.0 

2-coordinate of the surface data point. 

31-60 

DATA(4-6) FFSURF 
3F10.0 

(not used). 

Surface 

Data Property Card 

(6F10.0) 

1-10 

DATA(7) 

FFSURF 

F10.0 

Surface Mach number. 

11-20 

DATA(8) 

FFSURF 

F10.0 

X-direction cosine component of the surface 
velocity vector. 

21-30 

DATA(9) 

FFSURF 

F10.0 

Y-direction cosine component of the surface 
velocity vector. 

31-40 

DATA(10) 

FFSURF 

F10.0 

Z-direction cosine component of the surface 
velocity vector. 

41-50 

DATA(ll) 

FFSURF 

F10.0 

P/P . 

CIO 


51-60 DATA(12) FFSURF T/T . 






Surface Data Panel Selection Cards (1015) 

These cards are used when surface data are to be transferred from unit 4 to 
unit 10 by the routine FFSURF. The cards are input only if KTRNSF = 0 on 
the Flow Field Control Card. In subsequent streamline calculations the 
surface flow properties of each Panel's element centroids will be fit with 
the surface spline for interpolation purposes. 


Column 

Code 

Routine 

Format 

1-5 

IPANEL(I+1) 

FFSURF 

15 

6-10 

I PANEL(1+2) 

FFSURF 

15 

il-15 

IPANEL(1+3) 

FFSURF 

15 

16-20 

IPANEL(1+4) 

FFSURF 

15 

21-25 

I PANEL(I+5) 

FFSURF 

15 

26-30 

IPANEL(I+6) 

FFSURF 

15 

31-35 

I PANEL(I+7) 

FFSURF 

15 

36-40 

I PANEL(1+8) 

FFSURF 

15 

41-45 

I PANEL(1+9) 

FFSURF 

15 

46-50 

IPANEL(I+10)FFSURF 


Explanation 


Panel numbers of the surface data on unit 4 
to be transferred to the flow field storage 
unit 10. Panel numbers correspond to the 
order in which thn Panels were read by the 
geometry routines. 


NOTE: 


I =< 


0 

for 

1 < 

NSREG < 

10 

0,10 

for 

10 < 

NSREG < 

20 

0,10,20 

for 

20 < 

NSREG 1 

30 

0,10,20,30 

for 

30 < 

NSREG < 

40 

0,10,20,30,40 

for 

40 < 

NSREG £ 

50 


Example: for NSREG = 13, only 2 panel Selec¬ 
tion Cards are required. 
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2• SURFA CE STREAMLINE OPTION 


The new streamline method in the Mark IV program requires a different 
data format than that described on pp. 92-98 in Reference 1. The infor¬ 
mation presented here is intended to completely replace the old Surface 
Streamline option. 

The Surface Streamline option is called from the AERO executive pro¬ 
gram as described previously. The streamline method requires that the 
surface velocity direction cosines at all the Element centroids be de¬ 
fined a priori. There f ore, surface property information (M^, p/p^, T/T^, 
and the surface velocity direction cosines) for the entire vehicle must 
be available on unit 10 prior to accessing the Surface Streamline option. 
If the Inviscid Pressure option is used to generate surface data, the 
direction cosines of the surface velocity at each Element centroid are 
calculated using Newtonian theory (see Section III - Theory). If the data 
are hand-loaded using the Surface Data Transfer option, "exact" stream¬ 
lines may be traced if the surface velocity direction cosines are exact. 
The hand-loaded data points do not necessarily have to correspond with 
the Element centroids. 

One of the major problems associated with tracing streamlines over 
arbitrary Mark IV geometries is that of interpolating between the Element 
centroids for the surface flow properties. Since streamlines generally 
do not pass through the Element centroids, a means must be available for 
estimating the inviscid flow properties at any point on the surface. 

The Mark IV program employs the surface spline technique presented by 
Harder and Desmaris in Reference 3. The method generally cannot be 
applied to the entire vehicle surface since it is required that the quan¬ 
tity to be fit (e.g. p/p^.) be well-behaved within the boundaries of the 
spline. Therefore, the surface flow field must be divided into separate 
regions in which the surface properties do not exhibit rapid changes 
within each region. In the new streamline method, the surface properties 
corresponding to each geometry Panel are surface fit. For a geometry 
consisting of N Panels, a total of N independent spline fits are required, 
for each flow quantity. The six flow quantities that are fit include 
p/p , T/T , M, and the three components of the unit surface velocity. 

no oo 
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As defined in the Surface Data Transfer option, the flow properties 
associated with a Panel are termed a subregion. The use of Panel bound¬ 
aries for sub-dividing the surface flow field into well-behaved sub- 
regions is a logical approach since the inviscid properties, calculated 
using the local slope methods in the Mark IV program, do not exhibit 
radical behavior over a Panel unless the Panel itself contains rapid 
changes in character. Although it is difficult to concisely define a 
"well-behaved" flow property or geometry Panel, two rules of thumb may 
be used when preparing the geometry. First, a Panel should not contain 
rapid changes in curvature. For example, the longitudinal curvature of a 
spherically-blunted cone changes discontinuously from a nonzero value on 
the nosecap to a zero value at the sphere-cone juncture. The geometry 
should therefore be described by two Panels - the nosecap and conical 
frustum. Secondly, each Panel should contain only one geometry Section. 

Although the Mark IV geometry package permits the user to group several 
Sections into one Panel, no need exists for more than one Section within 
a Panel unless either the general shape or the boundaries of a surface 
change discontinuously. The curvature of the flat surface illustrated 
in Figure 3 is certainly well-behaved, but the slope of the boundaries 
changes discontinuously. The surface spline as used in the new stream¬ 
line method requires that each subregion, or Panel, consist of three or 
four boundaries whose slopes are continuous. Furthermore, since the 
spline is two-dimensional, each Panel must contain at least three Ele¬ 
ments whose centroids are not colinear (i.e., not a function of one 
spatial variable). 

In the new streamline method two options are available for initializing 
the streamline computation. First, the user may specify the points on the 
vehicle surface from which the streamlines are to be traced. Such points 
are termed starting points which may or may not correspond to true 
streamline origins. Since the user usually does not know the locations 
of the streamline origins a priori, an option has been included which 
allows the streamline equations to be integrated against the flow direction 
until a true origin is located by the streamline method. Running lengths 
saved along the streamline are then automatically reordered. For any 
given point on the body, therefore, this approach permits the calculation 

i 
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c. Incorrect means of grouping 
Elements: 1 Panel, 2 Sections 


Figure 3. Restriction on Geometry Preparation 




of the surface flow history prior to the point specified. If the user 
indicates that the equations are to be integrated in the direction of 
flow, the starting point is assumed to coincide with a true origin. 

The second option included in the new streamline method for initial¬ 
izing the streamline calculation allows the user to request that starting 
points be automatically distributed by the program. If this option is 
selected, starting points are distributed along the aft-most boundaries 
of the geometry and the streamlines are traced forward, against the flow 
direction, until the appropriate origins are located. After all running 
lengths are reordered, the method has the capability to determine 
whether the initial streamlines are sufficiently distributed for subse¬ 
quent skin friction calculations and integration. If the streamlines 
are not sufficiently distributed, new starting points are strategically 
distributed over the geometry surface and additional streamlines are 
computed. (Although the Viscous Methods option in the modified Mark IV 
program is capable of applying boundary layer methods along the surface 
streamlines, the present version of the code does not contain a method 
for assessing the contribution of viscous shear along the streamlines to 
total vehicle forces and moments. Instead, the user must resort to the 
Mark III Skin Friction Method included in the Viscous Methods option). 

Two criteria are used in the new streamline method for locating the 
origin of each streamline. The first criterion is simple: if the dot 
product of the freestream velocity and the local surface velocity is less 
than some small arbitrary value, an origin has been located. However, 
streamlines may originate from the leading edges of swept wings, nacelles, 
etc. in which the above criterion is useless. Therefore, as part of the 
input data preparation for the Surface Streamline option, the user must 
indicate whether each Panel contains a leading edge boundary, a trailing 
edge boundary, or all interior boundaries. A given Panel contains a 
leading edge boundary if no other Panels lie forward of and adjacent to 
the Panel. Trailing edge boundaries and interior boundaries are defined 
similarly. Provision has also been made for Panels that contain both 
leading and trailing edge boundaries. 
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In order to gain familiarity with the new streamline method, it is 
suggested that the user begin with a simple geometry such as a cone, a 
flat plate, or only a few Panels of the vehicle of interest. After the 
preliminary streamline calculation, other Panels may be easily added. 

The preparation of the input data for the Surface Streamline option 
is particularly simple. The data consists of four sets of cards: 

(1) Surface Property Access Card - indicates where the surface 

property data are to be found on unit 10. This card is 
read only once per entry into the Streamline option. 

(2) Streamline Data Card - contains general streamline informa¬ 

tion including flags that indicate which options are to 
be used. This card is also read only once per entry into 
the Streamline option. 

(3) Panel Description Cards - each card contains information 

regarding the general shape of a Panel. The number 
of these cards must equal the number of Panels on 
unit 4. The order in which the cards are input must 
coincide with the order in which the Panel information 
was stored on unit 4. 

(4) Streamline Starting Point Cards - indicate the points from 

which the streamline calculations are to begin. These 
cards follow the Panel Description Cards, and are used 
only if the starting points are to be input by the user, 
as specified on the Streamline Data Card. 

A detailed description of the formats for each set of cards is pre¬ 
sented in the following pages. 
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Column Code 


Routine 

Format 


Explanation 


1-2 

NOSET 

STREAM 

12 

Data set number where surface properties will 
be found on unit 10. 

3-4 

IABSET 

STREAM 

12 

a-8 set number where surface properties will 
be found on unit 10. 

5-6 

IREGON(l) STREAM 

12 

Flow region number where surface properties 
will be found on unit 10 if data were trans¬ 
ferred from unit 4. 


=1 If data were transferred from unit 4 to 
unit 10. 

= 0 If not. 

7-8 IREG0N(2) STREAM Flow region number where surface properties 

12 will be found on unit 10 if data were hand 

loaded. 

=2 If data were hand load. 

=0 If not. 

NOTE: The surface flow property information 
of the two regions taken collectively must 
uniquely define the surface properties of 
the entire vehicle. 


Streamline Data Card (512, 10X, 2F10.0) 

1-2 NSTR STREAM Total number of streamlines to be traced if 

12 ISTART f 2. If ISTART = 2, NSTR is the number 

of streamlines per aft Panel (KBNDRY = 2 or 
3 on the Panel Description card pp. 20) to be traced. 
No maximum allowable value is established for this 
parameter. 

3-4 IPRINT STREAM Print coordinates and corresponding flow pro¬ 

perties of each streamline. 

= 0 Do not print. 

* 1 Print. 
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Streamline Data Card (Continued) 


Column 

Code 

Routine 

Format 


Explanation 

5-6 

I STORE 

STREAM 

12 

Save information along each streamline on units 50 
and 51 (see Section IV-New Local Storage Units). 




= 0 

Do not save. 




= 1 

Save. 

7-8 

I START 

STREAM 

12 

Streamline starting condition flag. 

= 0 Start streamline calculation at the 
centroid of the given Element number 
of the specified Panel number. (See 
Streamline Starting Point Cards, pp. 21). 




= 1 

Start streamline calculations at the 
given X, V, Z locations. (See 

Streamline Starting Point cards, pp. 21). 




= 2 

Appropriate starting points will be 
distributed by the code. 

9-10 

MORSTR 

STREAM 

12 

Additional streamlines flag. Used only if 

I START = 2. 




= 0 

No additional streamlines. 




= 1 

Locate appropriate additional starting 
points to ensure that a sufficient 
streamline distribution exists for 
viscous force computations and integration. 

21-30 

DIRECT 

STREAM 

F10.0 

Specifies the direction of the streamline 
integration. 




= 1 

Integrate in the direction of the flow. 




= -1 

Integrate against the flow direction. 

The streamlines will be traced until a 
true origin is reached. The running 
lengths of the streamlines will be re¬ 
ordered so that the maximum running 
lengths occur at the starting points and 
a zero running length corresponds to the 
true origins. If ISTART = 2, DIRECT is 
automatically set to -1. 
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Streamline Data Card (Continued) 


Column 

Code 

Routine 

Format 

31-40 

XNOSE 

STREAM 

F10.0 


Explanation 


The axial location, X, of the forward-most 
point on the vehicle. 


Panel Description Cards (212) 

The number of these cards must equal the total number of subregions (geom¬ 
etry Panels) in regions IREGON(l) and IREG0N(2). The order in which these 
cards are read must correspond to the order in which the Panels were read 
by the geometry routines, e.g. the 3rd Panel Description Card must apply 
to Panel number 3. 

Column Code Routine Explanation 

_ __ Format _____ 

1-2 KSHAPE(I) STREAM Surface spline flag. 

12 

= 0 The Ith Panel is not a body of revo- 
1ution. 


= 1 The Ith Panel is best described in 
cylindrical coordinates. The spline 
will use the functional form R = f(A,q). 
KSHAPE(I) should be set to 1 only if the 
circumferential angles of all cross- 
sections are 180°. If the Panel is a 
true body of revolution (360°), the Panel 
must be divided into two Panels. 


3-4 KBNDRY(I) STREAM Panel boundary flag. 

12 

= 0 The Ith Panel is an interior Panel. Adja¬ 
cent Panels lie on all four sides (an 
adjacent Panel may include the plane of 
symmetry). 

= 1 The Ith Panel is a leading edge Panel. 

No Panels are adjacent to the forward 
most boundary of this Panel. 

= 2 The Ith Panel is a trailing edge Panel. 

No Panels are adjacent to the aft most 
boundary of this Panel. 
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Panel Description Cards (Continued) 


Column Code Routine 
Format 


Explanation 


3-4 KBNDRY(I) STREAM 

12 


= 3 The Ith Panel contains both leading and 
trailing edges. 


Streamline Starting Point Cards (213, 4X, 3F10.0) 

These cards are used only if ISTART = 0 or 1, see pp. 19. The number of these cards 
must equal NSTR, the number of streamlines to be traced. 


1-3 

LPANEL 

STREAM 

13 

The Panel number on unit 4 for the start of 
the streamline. Used if ISTART = 0 or 1. 

4-6 

L 

STREAM 

13 

Element number in Panel LPANEL for the start of 
the streamline. Used only if ISTART = 0. 

11-20 

XS 

STREAM 

F10.0 

X-coordinate of the streamline starting point. 
Used only if ISTART = 1. 

21-30 

YS 

STREAM 

F 10.0 

Y-coordinate of the streamline starting point. 
Used only if ISTART = 1. 

31-40 

ZS 

STREAM 

F10.0 

Z-coordinate of the streamline starting point. 
Used only if ISTART = 1. 
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3. VISCOUS PROGRAM OPTION 


The Viscous Program option, called from the AERO executive program 
as shown in Figure 1, is used for all viscous calculations. Two sub¬ 
options are available in the Viscous option for estimating the effects of 
the boundary layer. One sub-option is used when detailed boundary layer 
information is required at specific points on the body, and another is 
used to estimate the contribution of skin friction to total vehicle forces 
and moments. The following user-oriented descriptions of both sub¬ 
options are intended to completely replace the documentation given on 
pp. 114-128 in Reference 1, Volume I. 

The first sub-option applies boundary layer methods along inviscid 
surface streamlines generated previously by the Surface Streamline option. 
Selection of this sub-option requires that the streamline data be avail¬ 
able on unit 50 (see ISTORE on the Streamline Data Card described in the 
Surface Streamline option, pp. 19). Although Mark III Skin Friction methods(flat 
plate methods) may be used along the streamlines, the primary purpose of the 
current sub-option is to provide the user with a means for estimating de¬ 
tailed boundary layer properties along the streamlines. Integral boundary 
layer methods (References 4&5) are employed to provide such information 
as momentum thickness, displacement thickness, and velocity profiles as 
well as skin friction and heat transfer. 

If the user requires the contribution of skin friction to total 
vehicle forces and moments, the second sub-option must be used. This sub¬ 
option is known as the Mark III Skin Friction method, and employs flat 
plate methods to each of the Elements in the geometry and therefore does 
not provide variable properties along the streamlines. Since running 
lengths to eachof the Elements must be input, the user normally prepares 
a simplified geometry model. Prior to accessing the Mark III Skin Friction 
method it is necessary that inviscid surface property data, generated by 
the Inviscid Pressures option of the AERO program, be available. 

Modifications to the Viscous Program option were necessary to ensure 
that the integral boundary layer methods, which must be applied along 
inviscid surface streamlines, are compatible with the new streamline 
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method. Corrections were also made to the equations used in the integral 
methods as described in Section III - Theory. However, only minor changes 
were made to the overall input data procedure for the Viscous Program 
option. Although modifications were made to the initial data cards used 
by both sub-options, no changes were made to the format of the Mark III 
Skin Friction Element Data Cards. 


Viscous Method Card (12, 8X, 15A4) 

Column Code Routine Explanation 

__ Format _____ 

-2 ISFMTH VISCUS Viscous method flag. 

12 

= 0 Apply viscous methods along inviscid sur¬ 
face streamlines. Wall temperature, if 
it is not input, will be calculated by the 
Mark III skin friction methods. The 
Boundary Layer Method Control Card, 
pp. 24, will be expected next. The 
Mark III Skin Friction cards will 
not be input. 

= 1 Calculate skin friction coefficients 
using the Mark III program methods. 

The Mark III Skin Friction Coefficients 
Basic Flag Card, pp. 30, will be ex¬ 
pected 'next. 


ii'/'J TITLE VISCUS Title to be printed on the skin friction output 
15A4 pages. 
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INPUT DATA FOR VISCOUS CALCULATIONS ALONG STREAMLINES 

Viscous calculations along streamlines may be made only when ISFMTH = 0 
on the Viscous Method Card. If ISFMTH= 1, skip to the Mark III Skin Friction 
Coefficient Basic Flap Card, pp. 30. One of two methods may be selected for 
computing viscous effects along the surface streamlines. The more commonly 
used option employs integral boundary layer methods which are applicable 
to arbitrary pressure gradients. The second option consists of the same 
methods used by the Mark III Skin Friction program. Such methods are 
strictly applicable only to zero pressure gradient flows. 


Boundary Layer Method Control Card (12, 211, 12, II, 3X, 4F10.0) 


Column 

Code 

Routine 

Format 

Explanation 

1-2 

ISTRML 

INTEG 

12 

Streamline number for this set of viscous cal¬ 
culations. N Boundary Layer Method Control 
cards for N streamlines. No maximum value is 
assigned to this parameter. 

3 

LASTSL 

INTEG 

11 

Last streamline flag. 

= 0 This is not the last Boundary Layer 

Method Control Card. If ISFM = 0, 
another Boundary Layer Method Control 

Card is expected after the Integral 

Method Flag Card, pp. 27-29. If ISFM = 1, 
another Boundary Layer Method Control 

Card will be expected next. 




= 1 This is the last Boundary Layer Method 
Control Card. The program will return 
to the AERO routine after viscous cal¬ 
culations are made along this streamline. 

4 

ISFM 

INTEG 

11 

Boundary Layer Method Card for this stream¬ 
line. 


= 0 Use integral boundary layer method. 

= 1 Use one of the Mark III Skin Friction 
methods. Method will be selected by 
using the IWT flag below. 
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Boundary Layer Method Control Card (Continued) 


Column Code 


5-6 IWT 


Routine Explanation 

Format 


INTEG Wall temperature method flag. This flag 

12 controls the selection of the method to be 

used in calculating the wall temperature in 
routine TEMP. This flag is used for both 
ISFM = 0 and = 1. When ISFM = 1 it also 
controls the skin friction coefficient cal¬ 
culation procedure selection. In the dis¬ 
cussions below the methods to be used for 
laminar and turbulent flow are separated by 
a slash (i.e., Laminar/Turbulent). 

= 0 Use Reference Temperature/Spalding- 
Chi methods to calculate temperature. 

= 1 Use adiabatic wall temperature and 
Reference Temperature/Spalding-Chi 
methods. 

= 2 Use input wall temperature and Ref¬ 
erence Temperature/Spalding Chi meth¬ 
ods. Wall temperature is input in 
CC 11-20 and CC 21-30. 

= 3 Use Reference Enthalpy/Spalding-Chi 
(with enthalpy ratios) methods. 

= 4 Use adiabatic wall temperature and 

Reference Enthalpy/Spalding-Chi (with 
enthalpy ratios) methods. 

= 5 Use input wall temperature and Ref¬ 
erence Enthalpy/Spalding-Chi (with 
enthalpy ratios) methods. Wall 
temperature is input in CC 11-20 and 
CC 21-30. 

= 6 Use Reference Temperature/Reference 
Temperature methods. 

= 7 Use input wall temperature and Ref¬ 
erence Temperature/Reference Temperature 
methods. Wall temperature is input in 
in CC 11-20 and CC 21-30. 

= 8 Use Reference Enthalpy/Reference 
Enthalpy methods. 
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Boundary Layer Method Control Car d (Continued) 


Column Code Routine 
Format 


Explanation 


5-6 

IWT 

(Contin¬ 

ued) 

= 9 Use input wall temperature and Ref¬ 
erence Enthalpy/Reference Enthalpy 
methods. Wall temperature is input 
in CC 11-20 and CC 21-30. 

7 

IPRINT 

INTEG 

11 

Iteration and local skin friction print flag 
for use in routine TEMP. 




=0 Do not print. 




= 1 Print iteration results for wall tem¬ 
perature and the final local skin- 
friction data in routine TEMP. 




= 2 Print the final local skin-friction data 
from routine TEMP but do not print the 
iteration results. This is the recom¬ 
mended option for most applications. 

11-20 

SURFI6 

INTEG 

F10.0 

Input wall temperature for laminar calcula¬ 
tions, °R. This input is used when IWT = 

2, 5, 7, or 9. 

21-30 

SURFI7 

INTEG 

F10.0 

Input wall temperature for turbulent cal¬ 
culations, °R. This input is used when 

IWT = 2, 5, 7, or 9. 

31-40 

RETRAN* 

INTEG 

F10.0 

Transition Reynolds number X 10’^. Turbu¬ 
lent boundary layer calculations begin when 
the local edge Reynolds number attains this 
value. If ISFM = 0 and RETRAN = 0, the 
initial displacement thickness (DTURB) and 
the initial momentum thickness (TTURB), 
specified on the Integral Method Flag Card 
below, must both be non-zero. 

41-50 

STRAN* 

INTEG 

Transition running length, ft. Turbulent 
boundary layer calculations begin when the 
running length along the streamline attains 
this value. If ISFM = 0 and STRAN = 0., 
the initial displacement thickness (DTURB) 
and the initial momentum thickness (TTURB), 
specified on the Integral Method Flag 

Card below, must both be non-zero. 


transition occurs when either RETRAN or STRAN is exceeded. 





This card is input only when ISFM = 0. 


Column 

Code 

Routine 

Format 

Explanation 

1-2 

NVP 

INTEG 

12 

Number of points desired in the velocity pro¬ 
file at each station. (Usually input = 20). 

No maximum value is assigned to this parameter. 

3-4 

KSMTH 

INTEG 

12 

Number of times distribution of surface veloc¬ 
ity is to be smoothed prior to computation of 
surface gradients ( = 0,1,2,3, etc.). No maximum 
value is assigned to this parameter. 

5-6 

KSPLN 

INTEG 

12 

Integer indicating manner in which surface 
gradients are to be calculated. 

= 0 Weighted-difference technique (preferred). 

= 1 Spline curve- f it technique. 

7-8 

KLE 

INTEG 

12 

Flag indicating the type of initial condi¬ 
tion existing at the first streamline point. 

= 0 Stagnation point or initial values 
given. 

= 1 Sharp leading edge. 

9-iO 

KATCH 

INTEG 

12 

Flag indicating whether laminar boundary 
layer separation (if encountered) should 
reattach as a turbulent boundary layer. 

= 0 Stop calculations if separation 
encountered. 

= 1 Reattach. The input variable CTHET 
(below) must be non-zero. 

11-12 

KPRE 

INTEG 

12 

Preliminary calculation print flag. 

= 0 Output suppressed. 

= 1 Output printed. 

13-14 

KGRAD 

INTEG 

12 

Print flag for surface velocity and Mach 
number gradients. 


= 0 Output suppressed. 


= 1 Output printed. 
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Integral 

Method 

Flag Card 

(Continued) 

Column 

Code 

Routine 

Format 

Explanation 

15-16 

KSDE 

INTEG 

12 

Flag for printing direct results of the 
integration of the laminar and turbulent 


boundary layer equations. (Usually input = 
0 ). 

= 0 Output suppressed. 

= 1 Output printed. 

17-18 KLAM INTEG Flag for printing results of laminar insta- 

12 bility and transition calculations. 

= 0 Output suppressed. 

= 1 Output printed. 

19-20 KMAIN INTEG Flag for printing principal boundary layer 

12 information (C^, Nu, 6*, 0, etc.). Usually 

input = 1. 

= 0 Output suppressed. 

= 1 Output printed. 

21-22 KPROF INTEG Flag for printing velocity profiles. 

12 

= 0 Output suppressed. 

= 1 Output printed. 

31-40 CTHET INTEG Ratio of momentum thickness after reattach- 

F10.0 ment to momentum thickness at laminar 
separation. This parameter used when 
KATCH = 1 if separation occurs. 

41-50 DLAM INTEG Initial displacement thickness (ft.), if 

F10.0 any, of the laminar boundary layer. If 

DLAM is zero, initial laminar displacement 
thickness will be calculated by the program 
according to the value of the KLE flag 
specified above. 

51-60 TLAM INTEG Initial momentum thickness (ft.), if any, 

F10.0 of the laminar boundary layer. If TLAM is 
zero, initial laminar momentum thickness 
will be calculated by the program according 
to the value of the KLE flag specified above. 
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Explanation 


i 


I n teg ra 1 Method Flag Card (Cont i n ued ) 

Column Code Routine 
Format 


61-70 DTURB INTEG 
F10.0 


71-80 TTURB INTEG 
F10.0 


Initial displacement thickness (ft.), if 
any, of the turbulent boundary layer. If 
the boundary layer is turbulent at the 
first streamline point (RETRAN = 0. or 
STRAN =0.), DTURB must be non-zero. If 
transition occurs downstream of the first 
streamline point, DTURB may be set to zero 
and the laminar value at the transition 
point is used. 

Initial momentum thickness (ft.) if any, of 
the turbulent boundary layer. If the bound¬ 
ary layer is turbulent at the first stream¬ 
line point, TTURB must be non-zero. If 
transition occurs downstream of the first 
streamline point, TTURB may be set to zero 
and the laminar value at the transition 
point is used. 


NOTE: If LASTSL = 0, another Boundary Layer Method Control Card will be 
expected after the above card. If LASTSL = 1, the program will return to 
AERO. 
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MARK III SKIN FRICTION METHOD 


Mark III Skin Friction Basic Flag Card (12, 211) 


Column Code 


NCOMP 


IFSAVE 


Routine 

Format 

VISCUS 

12 


VISCUS 

II 


Explanation 


Total number of vehicle Components to be 
analyzed. Each Component may consist of 
one or more vehicle panels. The grouping 
of Panels to form Components is controlled 
by the Geometry Data Source Card below. 

Force data save flag. 

= 0 Set up a new force data save file 
(unit 9). Save skin friction force 
data for future summation. 


1 Save skin friction force data on 
unit 9 for future summation. Use 
old unit 9 file and just add the new 
force data onto the file. 

2 Do not place the force data on the 
force data file unit. 


IPRINT 


VISCUS 

II 


Skin friction print flag. 

= 0 Do not print. 

= 1 Print detailed skin friction inter¬ 
mediate results. 


Geometry Data Source Card (2012, IX, 13) 


IPANL(l) VISCUS 
2012 


39-40 IPANL(20) 
41-44 NS 


VISCUS 

13 


The identification numbers for all of the 
Panels on the Quadrilateral Element Storage 
unit (4) that are to be grouped to form 
this vehicle component. 


Number of skin friction elements to be 
analyzed. This number must be equal to 
the number of elements on the Quadrilateral 
Element save unit (4) for this vehicle 
Component and must not be greater than 100. 
The number of Skin Friction Element Data 
Cards must be = NS. This input is used 
for the Mark III skin friction option only. 
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Mark III Skin Friction Element Data Cards 
(12, 811, 2F9.0, 3F6.0, 2F6.0, F4.0, 8X, 12) 


One Skin Friction Element Data Card must be loaded for each element stored 
on the Quadrilateral Element Storage unit (4) for each vehicle Component. 
The format of these cards is exactly the same as the Type 11 cards used 
on the Mark III program (Mode 1 skin friction method). However, some of 
the parameters on the old Type 11 card are not actually used by this new 


version 

of the 

program. 


Column 

Code 

Routine 

Format 

Explanation 

1-2 

IS(I,1) 

SKINFR 

12 

Skin friction element number. 

3 

IS(I,2) 

M3SF 

11 

Viscous-Inviscid interaction effect flag 

= 0 Use tangent-wedge in interaction 
correction. 




= 1 Use tangent-cone in interaction 
correction. 

4 

IS(I,3) 

SKINFR 

11 

Calculate induced pressures due to bound 
ary layer displacement effects. Skin 
friction is not calculated. 




= 0 No 




= 1 Yes 

5 

IS (1,4) 

SKINFR 

11 

Skin-friction summation flag. 


= 0 Use turbulent skin friction data in 
calculating forces. (Note: The pro¬ 
gram will make a switch to laminar 
summation at very low Reynolds number, 
where turbulent results are not 
meaningful). 

= 1 Use laminar skin friction data in 
calculating forces. 

6 IS(1,5) SKINFR (Not used in this program). 

II 
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Mark III Skin Friction Element Data Cards (Continued) 


Column 


7 


Code 


IS{I,6) 

laminar and turbulent skin-friction results 
The result to be added to the pressure cal¬ 
culations is indicated by the flag in CC 5. 
In the discussions below the methods to be 
used for laminar and turbulent flow are 
separated by a slash (i.e., Laminar/Turbu¬ 
lent). (Integer) 

= 0 Calculate wall temperature and skin 
friction using Reference Temperature/ 
Spalding-Chi methods. 

= 1 Use adiabatic wall temperature and 
Reference Temperature/Spalding-Chi 
methods. 

= 2 Use input wall temperature and Refer¬ 
ence Temperature/Spalding-Chi methods. 

T input in CC 47-52 and 53-58. 
w 

= 3 Calculate wall temperature and skin 
friction using Reference Enthalpy/ 
Spalding-Chi (with enthalpy ratios) 
methods. 

= 4 Use adiabatic wall temperature and 
Reference Enthalpy/Spalding-Chi 
(with enthalpy ratios) methods. 

= 5 Use input wall temperature and Refer¬ 
ence Enthalpy/Spalding-Chi (with 
enthalpy ratios) methods. T input 
in CC 47-52 and 53-58. w 

= 6 Calculate wall temperature and skin 
friction using Reference Temperature/ 
Reference Temperature methods. 

= 7 Use input wall temperature and Refer¬ 
ence Temperature/Reference Temperature 
methods. T input in CC 47-52 and 
53-58. W 


Routine 

Format 

SKINFR 

II 


Explanation 


Wall-temperature and skin-friction method 
Flag. The program always calculates both 



Mark III Skin Friction Element Data Cards (Continued) 

Column Code Routine Element 

_ _ Format ___ 

7 IS(1,6) (Continued) = 8 Calculate wall temperature and skin 

friction using Reference Enthalpy/ 
Reference Enthalpy methods. 

- 9 Use input wall tempera-d Refer¬ 
ence Enthalpy/Reference fit'' lpy 
methods. T input in CC and 

53-58. w 

8 IS(I»7) SKINFR Flag to control printing of skin-friction 

II data for each skin-friction surface element. 

=0 Do not print. 

= 1 Print skin-friction data. This is 
recommended option for most applica¬ 
tions. 

9 IS(I,8) SKINFR Print flag for flow characteristics before 

II and after the shock or expansion. 

= 0 Do not print. 

= 1 Print flow characteristics. 

10 IS{1,9) SKINFR Iteration and local skin friction print 

II flag. 

= 0 Do not print. 

= 1 Print iteration results for wall 
temperature and the final local 
skin-friction data. 

= 2 Print the final local skin-friction 
data but not the iteration results. 

This is the recommended option for 
most applications. 

11-19 SURF(1,1) SKINFR Skin friction element surface wetted area in 

F9.0 same units as S re f. If input as 0.0 then 

the program will use the surface area as cal¬ 
culated from the input geometry unit for each 
element. The input wetted area must corre¬ 
spond to the input skin-friction geometry 
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Mark III Skin Friction Element Data Cards (Continued) 

Column Code Routine Explanation 

Format 


11-19 SURF(I,1)(Continued) (i.e., if the Symmetry flag is 0, left side 

of the vehicle input, then the input wetted 
area should be only for the left side). 

The four input quantities in CC 20 through 46 furnish to the program the 
planform shape of the skin-friction surface being analyzed ("Surface-of- 
Interest"), and the shape of the initial-surface (to account for the fact 
that the flow has traversed some other part of the shape before reaching 
the surface of interest). This information is not obtained from the 
input skin-friction geometry data input on the Type 3 cards. The input 
skin-friction geometry data are used only to establish the position and 
orientation of the centroid and the area of each skin-friction surface. 
The diagram below illustrates the input parameters required on the Skin 
Friction Element Data Cards. 



Initial -Surface 


Surface-of-Interest 







Hark III Skin Friction Element Data Cards (Continued) 


Column Code Routine 
_ Format 


Explanation 


35-40 SURF(I,4) SKINFR 
F6.0 


41-46 SURF(1,5) SKINFR 
F6.0 


The taper ratio of the initial-surface 
Ul/Li). This taper ratio is defined as 
the ratio of the shortest chord length to 
the longest chord length. If both the 
initial-surface longest-length and the 
longest length of the surface-of-interest 
are on the same edge of the shape, then 
the taper ratio of the initial-surface is 
input as a positive number. If these lengths 
are on opposite sides of the shape such as 
in the diagram on the previous page then the 
initial surface taper ratio is input as a 
negative number. With these ground rules 
the absolute value of the taper ratio will 
never be greater than 1.0. 

The taper ratio of the surface-of-interest 
(& 2 /L- 2 )' This taper ratio is defined as 
the ratio of the shortest chord length. 

This taper ratio is always positive and 
never greater than 1.0. 


47-52 

SURF(1,6) 

SKINFR 

F6.0 

Input wall temperature for laminar calcu¬ 
lations, °R. This input is used when CC 7 = 
2, 5, 7, or 9. 

53-58 

SURF(I,7) 

SKINFR 

F6.0 

Input wall temperature for turbulent calcu¬ 
lations, °R._ This input is used when CC 7 = 
2, 5, 7, or 9. 

59-62 

SURF(I,8) 

SKINFR 

F4.0 

(Not used in present program) 

71-72 

TYPE 

SKINFR 

12 

Card Type number. Not used in present 
program. 
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CASE I: OGIVE-CYLINDER GEOMETRY 
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CASE 2: SPHERICALLY BLUNTED BICONE 
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EXPLANATION OF PRELIMINARY CALCULATIONS 
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tuRHOL Pul DIFFERtNUAL EOuATlONS - SOLUTION FOR F AND FOklil 
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F first dependent variable 'In turbulent boundary layer equations 

(related to momentum thickness) } See Section III - Theory 

FORMI second dependent variable, the Incompressible form factor 
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THET momentum thickness, 9 

DELTA boundary layer thickness, < 

FORM compressible form factor 

FORM! incompressible form factor 
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SECTION III 


THEORY 

Discussed in this section are the modifications made to the version 
of the Mark IV program released in November 1973 (Reference 1). The two 
primary changes made to the program involve surface streamline tracing and 
the application of integral boundary layer methods along the surface 
streamlines. 

The original version of the Mark IV program contains both a Newtonian 
streamline tracing method and a means for calculating detailed boundary 
layer properties along the streamlines. Due to the nature of the surface 
spline used in the streamline method to interpolate for the surface pro¬ 
perties, the surface flow field must be divided into regions in which the 
flow properties are "well-behaved," i.e., the properties do not exhibit 
sudden changes within a region. However, no means exists in the original 
version of the Mark IV program for tracing streamlines from region to region. 
Therefore, the arc lengths along the streamlines, necessary for subsequent 
boundary layer calculations, are correct only for those streamlines whose 
corresponding region contains a true origin. 

Other difficulties with the old streamline method are concerned with 
the location of the streamline origins. The user must supply the program 
with the origin for each streamline, but streamline origins are not usually 
known a priori. Furthermore, if the skin friction distribution is to be 
integrated to obtain the contribution of viscous shear to vehicle forces 
and moments, the starting points for the streamlines must be strategically 
placed to ensure that the resultant streamlines are sufficiently distributed 
over the geometry surface. Since streamline divergence is particularly 
evident on three-dimensional bodies, small changes in the placement of 
starting points normally result in large variations of the final streamline 
distributions, and the user cannot be expected to make a proper selection 
of the starting points. (Starting points used in this context are not 
necessarily true origins. Since the surface velocity field contains 
singularities at stagnation points, streamline calculations often begin 
some small distance e from the true origin.) 
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The old streamline method has been replaced by one capable of tracing 
continuous surface streamlines. In addition, streamline origins may be 
located automatically by the program. However, an option has been provided 
to allow the user to specify starting points if so desired. Logic has 
also been included which identifies those surfaces of the geometry that 
are not sufficiently covered by the initial streamlines. If more stream¬ 
lines are required, additional starting points are strategically placed 
and more streamlines are traced. 

The primary purpose of tracing surface streamlines is to provide paths 
for subsequent boundary layer calculations. The particular methods used 
in the Mark IV code are integral methods (References 4 & 5) originally 
coded by McNally (Reference 6) in a project unrelated to the Mark IV 
program. Since the integral methods rely on streamline informatin, the 
replacement of the old streamline method required minor modifications to 
the FORTRAN coding of the boundary layer methods. 

However, the results of some simple flat plate calculations revealed 
several discrepancies related to modifications made to McNally's FORTRAN 
program by the authors of the Mark IV program. As mentioned in Volume II 
of Reference 1, modifications to the McNally program were necessary to 
remove the assumption of isentropic flow used implicitly in the coding of 
the integral methods. However, the isentropic assumption was not properly 
removed, and steps were taken in the current effort to correct the problems. 
In addition, the integral equation used in the laminar boundary layer 
method (equation 32 of Reference 4) was replaced by the original form of the 
equation, a nonlinear, ordinary differential equation (equation 27 of Refer¬ 
ence 4). Given in the "Integral Boundary Layer Methods" sub-section of 
this report are brief descriptions of the integral methods used in the 
Mark IV program for the laminar and turbulent boundary layers. Changes 
made to the methods and to the FORTRAN coding of the methods are also 
discussed. 
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1. Streamline Tracing 

Although several methods are available for tracing inviscid surface 
streamlines {References 8 & 9), the one approach which is consistent with 
the engineering design methods of the Mark IV program is the Newtonian 
method, also known as the Steepest Descent method. The Newtonian model 
assumes that a stream of particles impinging on a surface retains its 
tangential component of momentum. Thus, the velocity at any point on the 
surface is assumed to lie in the plane formed by the freestream velocity 
and the local outward normal, 

V s = n x (V^x n) (1) 

where n is the unit outward normal and V is the freestream velocity. 

Since the surface velocity is simply the time rate of change of the position 
vector, r, the Newtonian streamline equation may be written 

| = n itVj n) (2) 

Using the definition of the magnitude of the surface velocity, V<. = 

dS/dt, where S is the streamline arc length, the vector streamline equation 
becomes 


dr 

dS 



(3) 


where V<. is the unit surface velocity calculated from Equation 1 . 

The vector Equation 3 is a system of three ordinary differential 
equations, but it is necessary to solve only two of the equations. Since 
the streamlines must lie on the geometry surface at all times, one coord¬ 
inate of the streamline position is related to the other two by the function 
describing the surface, X^ = f^, X^), where r = X i + X^j + X^k. 

One of the most difficult problems related to arbitrary-body streamline 
tracing is that of interpolating for the inviscid flow properties. Although 
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the unit surface velocity is of interest for Newtonian streamline calcu¬ 
lations, such edge conditions as static pressure and temperature are 
required for subsequent boundary layer calculations. In the Mark IV pro¬ 
gram the inviscid properties are calculated only at the Element centroids. 

A surface interpolation method must therefore be available to estimate 
the flow properties, including the unit surface velocity, between the 
centroids. 

The Mark IV program employs the surface interpolation method pre¬ 
sented by Harder and Desmaris in Reference 3. This method is not a 
parametric spline, i.e., an appropriate choice of the two independent 
variables is required to avoid multiple-values of the dependent variable. 

For example, attempts to surface fit a body of revolution in a Cartesian 
coordinate system would result in multiple-values of the dependent variable, 
regardless of whether the X-, Y-, or Z- coordinate is chosen as the depen¬ 
dent variable. It would be more appropriate to fit a body of revolution 
in cylindrical coordinates, R = R(A,t), where A is an axial distance and 
:■ is a meridonal angle. In the new streamline method, the surface spline 
for each Panel assumes one of three functional forms:Z =Z(X,Y), Y= Y(X,Z). or 
R = R(A,$). Parametric splines, however, do not require the appropriate 
selection of the two independent variables. 

The Harder and Desmaris surface spline, as used in the Mark IV pro¬ 
gram, has the advantage that the known coordinates used in generating the 
spline need not be ordered in a rectangular array, but has the disadvantage 
that it does not guarantee continuity of the fit between geometry panels. 

The spline must be applied to each geometry Panel, independently of all 
other Panels, and the quantity to be fit for each Panel (pressure, 
temperature, one appropriate geometry coordinate, etc.) must be well- 
behaved. Although it is difficult to succinctly define "well-behaved" in 
this context, experience has shown that the shape of the surface within the 
boundaries of a Panel must not contain rapid changes in curvature. Since 
local-slope pressure methods are used in the Mark IV program, the inviscid 
surface properties are generally well-behaved within a Panel's boundaries 
if the shape of the Panel is well-behaved. 
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The nonparametric form of the Mark IV surface spline and the lack of 
continuity of the surface fits between Panels seriously complicate the 
development of an arbitrary-body streamline method. To facilitate the 
streamline tracing a pseudo parametization of the surface spline is 
introduced. As shown in Figure 6, this parametization maps the 2-D domain 
associated with each Panel to a unit square, known as the u-v plane. A 
true parametric spline maps all three coordinates to the u-v plane. One 
advantage of working in the u-v plane is that the streamline integration 
is considerably simplified. After each integration step, a check must 
be made to determine if the streamline has crossed the boundaries of the 
current Panel. If the streamline integration is performed in the u-v 
plane, the check involves a simple test of the values of the independent 
variables, u and v: a streamline exceeds the boundaries of a given Panel 
if u or v is greater than 1 or less than 0. A similar check in Cartesian 
or cylindrical coordinates is far more complicated. Once a streamline 
exceeds the boundaries of a given Panel, the surface spline of the adjacent 
Panel must be used to continue the interpolation of the unit surface 
velocity required by the streamline equations. 

Once the appropriate functional form of a particular Panel is deter¬ 
mined, Y = Y(X,Z) or Z = Z(X,Y) or R = R(A,cj>), a bilinear mapping is used 
to transform the domain of the Panel's spline to the unit square, 

Xj(u,v) = X u (l-u)(l-v) + X 12 (l-u)v + X 13 uv + X 14 u(l-v) (4a) 

X 2 (u,v) = X 21 (l-u)(l-v) + X 22 (l-u)v + X 23 uv + X 24 u(l-v) (4b) 

where, as shown in Figure 7, X.j is the i^* 1 coordinate of the j^ corner 
point of the Panel's boundary, and X^ is the i^ coordinate of any point 
in the domain. Thus, the point (u,v) = (0,0) on the unit square corres¬ 
ponds to the point (X^,X 2 ) = ^ rom Equation 4 . Depending upon 

the appropriate functional form of the Panel, the coordinate pair (X^,X 2 ) 
represents (X,Z), (X,Y) or (A,4>). 

Once the appropriate domain of each Panel is mapped to the unit u-v 
square, the inviscid surface properties and the third geometry coordinate 
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of the centroids of each Panel are surface fit with u and v as the inde¬ 
pendent variables, i.e. p = p(u,v), T = T(u,v), V x = \l^( u,v), etc. One 
spline fit is required for each surface property of each Panel. The sur¬ 
face fits are calculated for all Panels prior to the streamline calculations, 
and the spline coefficients are stored on a random access unit for later 
use. 

If the streamlines are to be integrated in the u-v plane, the stream¬ 
line equations must reflect the change of coordinates. Choosing two of the 
components of the vector Equation 3 , the transformation to the u-v plane 
yields 
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and. 


I (V x ,V y ) , if (X 1§ X 2 ) = (X.Y) 

(V X ,V Z ) , if (X lf X 2 ) = (X,Z) 

(V A ,yR), if (x r x 2 ) = (A,d>) 

All partial derivatives may be evaluated from Equations 4a and b which 
describe the bilinear mapping. 

Given a starting point, the Panel associated with that point is 
identified, and the corresponding spline coefficients are retrieved from 
the random access unit. The starting point is transformed to the u-v 
plane, and the spline coefficients are used to interpolate for the sur¬ 
face properties at the starting point. A quartic Runge-Kutta integration 
of the streamline equations, 5a and b, is used to trace the streamlines 
across the u-v plane of the given Panel. At each step of the integration, 
u and v are mapped back to the Mark IV Cartesian coordinate system, and 
the (X,Y,Z) coordinates of the streamline are saved every N integration 
steps. 

If u or v is greater than 1 or less than 0 the streamline has exceeded 
the boundaries of the current Panel, and a search must be made for the 
adjacent Panel. Since the Panels need not be input by the user in any 
specific order, and since more than one Panel may be adjacent to one 
boundary of a Panel, a brute force search is performed to locate the proper 
adjacent Panel. Before searching for the adjacent Panel, the streamline 
is integrated an arbitrarily small distance e past the boundary of the 
last Panel. The Cartesian coordinates of the last streamline point are 
then mapped to the u-v planes of all remaining Panels. If the transformed 
point does not lie within the unit u-v boundaries of a given Panel, that 
Panel cannot be the desired adjacent Panel. However, if the transformed 
point does lie within a Panel's u-v boundaries, that Panel may or may not 
be the desired adjacent Panel. As shown in Figure 8, the domains of two 
Panels may overlap. To ascertain which Panel is the one of interest, the 
distance, d, between the streamline point and the surface of each of the two 
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Figure 8. Overlapping Spline Domains 
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Panels is computed. The Panel nearest the streamline point is assumed to 
be the desired adjacent Panel. Once the correct adjacent Panel is located 
the streamline integration continues. 

Integration of the streamline equations may be performed along or 
against the flow direction. Since the true origins of the streamlines are 
not known a priori, an option has been provided in the new streamline 
method which allows the starting points for the streamlines to be auto¬ 
matically distributed along the aft-most boundaries of the vehicle. The 
streamlines are then traced against the flow direction until the true 
origins are located. Since streamlines may originate from stagnation 
regions or from sharp leading edges, two criteria are used to determine 
when a streamline has reached its true origin: (1) the streamline 

/V A 

encounters a stagnation region, e, or (2) the streamline reaches 

a leading edge (a Panel boundary that is not adjacent to any upstream 
Panel). 

Starting points may also be input by the user. If the direction of 
integration is specified by the user to be against the surface flow, the 
procedure followed when starting points are program-specified is employed. 
If the streamlines are traced in the direction of flow, the starting points 
are assumed to be true origins, and the integration continues until a 
trailing edge or a staynation region is encountered. A trailing edge is 
any Panel boundary that is not adjacent to any downstream Panel. 

A technique has also been included in the new streamline method 
which strategically distributes additional starting points, based on the 
initial streamline distribution, to ensure that the streamlines are 
sufficiently distributed for future skin friction calculations and inte¬ 
gration. As the streamlines in the initial distribution are integrated 
across Panel boundaries, the intersections of the streamlines with the 
Panel boundaries are saved. In this manner, the streamlines' entrance 
and exit points to all Panels are known. A particularly simple algorithm 
may be developed if the points are saved in the u-v plane. As shown in 
Figure 9, straight lines are constructed in the u-v plane between the 
entrance and exit points of each Panel. If any of the areas bounded by 
the straight line segments and a Panel's boundaries is greater than some 
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Figure 9. Algorithm Used to Ensure Sufficient Streamline Distribution 
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arbitrary value (.30 used presently), appropriate additional starting points 
are distributed. 

The algorithm is applied to all Panels of the geometry starting with 
the aft-most Panels. In this manner the nunber of additional streamline 
computations is kept to a minimum since the additional streamlines from 
the aft Panels tend to fill-in the distribution over the upstream Panels. 

If an upstream Panel requires additional streamlines, the integrations 
are performed only in the forward direction, toward the streamline origins. 
Additional streamlines are never traced in the aft direction. 

Validation of the FORTRAN coding of the Newtonian streamline equations 
was straightforward. A cylindrical geometry at angle of attack was used 
in the validation since the outward normal at any point may be expressed 
analytically. Simple hand calculations of the unit surface velocity at 
points along the streamlines were in excellent agreement with the corres¬ 
ponding streamline directions. 

Comparisons of the Newtonian method (Steepest Descent) with more 
sophisticated inviscid streamline methods have been made by other investi¬ 
gators (Reference 9 is one among many). No such detailed comparisons 
were part of the present effort. 

Shown in Figure 5 is the streamline distribution over the X-24C 
forebody at 5 degrees angle of attack. The streamlines were generated 
using the program-specified starting point option in the new Mark IV 
streamline method. The plot was produced by the TEKtronix Pictures 
(TEKPIC) program (Reference 2), an interactive tektronix picture drawing 
code developed as part of the present effort. 
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2. Integral Boundary Layer Methods 

Given the inviscid condition' long the three-dimensional streamlines, 
boundary layer methods may be employed to predict such quantities as skin 
friction and heat transfer. The boundary layer methods chosen by the 
original authors of the Mark IV program include the method due to Cohen 
and Reshotko (Reference 4) for the laminar boundary layer, and the method 
presented by Sasman and Cresci (Reference 5) for the turbulent boundary 
layer. Both methods are integral approaches applicable to two-dimensional, 
compressible flows with arbitrary pressure gradients and heat transfer. 

The Schlichting-Ulrich method (Reference 10) is used to predict the point 
of neutral stability of the laminar boundary layer, and the distance be¬ 
tween the point of instability and the transition point is predicted using 
thn empirical curve presented by Granville (Reference 11). 

Each of the methods were incorporated into a FORTRAN computer code by 
McNally (Reference 6), and later modified for use in the Mark IV program. 
During the current effort, initial modifications were made to the peripheral 
coding (COMMON blocks, random access read/writes, etc.) to insure that the 
integral methods received the proper inviscid edge conditions from the new 
streamline method. However, the results of subsequent flat plate calcula¬ 
tion using the integral methods, as coded in the Mark IV program, were in 
disagreement with accepted simplified theories (Eckert's reference enthalpy 
and Van Driest II). As described below, the errors were due primarily to 
the fact that several of the variables appearing in the equations were 
based on freestream conditions rather than local conditions. In addition 
to correcting such errors, the current effort involved the replacement of 
Cohen and Reshotko's linearized equations with the nonlinear form of their 
equations. 

Only a brief description of the integral methods is presented here. 

The derivations of the laminar and turbulent boundary layer method may be 
found in References 4 and 5, respectively. Detailed descriptions of the 
various curve fits used in the program (e.g., thermal conductivity, 
k = k[T]), the integration schemes, and the FORTRAN coding of the methods 
are given in Reference 6. 
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a. Laminar Boundary Layer Method 

Cohen and Reshotko's approximate method is based on an approach first 
introduced by Thwaites. The standard Stewartson-Illingsworth transforma¬ 
tion is first applied to Prandtl's compressible boundary layer equations, 
and the transformed equations are then integrated with respect to the 
coordinate normal to the surface to yield a set of nonlinear, first order 
differential equations--the integral equations. These equations are then 
expressed in terms of three dimensionless parameters related to wall 
shear, surface heat transfer, and the pressure gradient. The momentum 
equation becomes. 


-U 

e 



where. 


2[n(H tr + 2) + i] 
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A similar procedure may be followed for the energy equation, but the 
following assumption negates the need for the energy equation. Analogous 
to Thwaites approach, universal functions were sought for the shear 
parameter, £, and the heat transfer parameter, r, in terms of the pressure 
gradient parameter, n, and the wall enthalpy function, S w = jp - 1. Such 
relationships were extracted from similar solutions (Falkner-§kan type 
flows) to the compressible, laminar boundary layer equations (Reference 12). 
Therefore, the fundamental equation to be solved is 


-u 4-U 

e dX U 


= N(n,S ) 
' w' 


where N(S w> n) is given from similar solutions. Given n along a streamline, 
the parameters ? and r may be determined from the curve fits 2 = «.(n,S w ) 
and r = r(n,S ) obtained from the similar solutions. Finally, the momentum 
thickness, the skin friction, and the heat transfer may be calculated from 
the definitions of n, Z, and r, respectively. 

Cohen and Reshotko simplified Equation 7 by noting that the function 
N(n,S w ) is approximately linear in n for fixed values of S w , 

N = A + Bn (8) 


Equation 7 may then be solved to yield, 
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or, in terms of physical quantities. 



where 

K = (3y - l)/(2y - 2). 

The latter equation was coded into McNally's FORTRAN boundary layer 
program which was later modified for use in the Mark IV program. In 
Volume II of Reference 1 it was stated that the isentropic assumption used 
implicitly throughout the boundary layer equations required a major modi¬ 
fication of McNally's coding. This modification primarily consisted of 
replacing (T /T ), wherever it appeared in the equations, by (p /p e )^ Y 
It is not clear why this thermodynamic relation was used. The isentropic 
edge assumption, mentioned by McNally in Reference 6, applies to the man¬ 
ner in which the edge conditions are calculated. Since McNally was con¬ 
cerned only with the boundary layers of shock-free flowfields, such 
quantities as the kinematic viscosity, v q = P Q /p 0 » appearing in the defi¬ 
nition of the pressure gradient parameter n, were based on freestream 
conditions. However, since the Mark IV program was designed for supersonic- 
hypersonic flows, the entropy changes instantaneously across the bow shock, 
and continues to vary along each surface streamline (edge temperature in 
the Mark IV code is calculated by assuming tangent-wedge conditions to 
exist locally). Therefore, by using local conditions instead of freestream 
conditions, the "isentropic assumption" is removed. 

The other modification made to the Mark IV integral methods involved 

the laminar boundary layer equation itself. As may be observed from the 

plot of N(n,S ) versus n with S as the parameter in Figure 4 of Reference 4, 
w w 

N is linear in n only for values of S very near zero, i.e., for approx- 

w 

imately adiabatic walls. To circumvent this problem, McNally allowed the 
constants A and B to vary with the pressure gradient parameter n. However, 
the mathematical validity of fixing coefficients in a differential equation, 
solving the equation, and subsequently allowing the coefficients to vary 
is questionable. Therefore, the linearized equation was replaced with 
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the original nonlinear relation. Equation 7 . Transforming this equation 
back to physical quantities there results 



N(n,S ) 
' w 



( 11 ) 


Equation 11 can easily be solved using the quartic Runge-Kutta 

numerical integration technique. The initial value of n depends on whether 

the streamline originates at a stagnation region or at a sharp leading 

edge. Curve fits of n_ p versus S obtained from similar solutions are used 

w 

for stagnation regions. If the streamline originates at a sharp leading 
edge, n = 0 . 


b. Turbulent Boundary Layer Method 

The method due to Sasman and Cresci employs the momentum integral and 

the moment-of-momentum integral equations for arbitrary pressure gradients. 

A Mager-Type transformation is then used to simplify the equations. Rather 

than solve the energy equation simultaneously with the momentum and moment- 

of-momentum equations, Sasman and Cresci assumed the Crocco relation to 

hold for flows with heat transfer and pressure gradient. In addition, the 

Ludwieg-Tiliman skin friction relation for incompressible turbulent flow 

(Reference 13 ) was used to relate the incompressible skin friction 

coefficient, C*, to the transformed adiabatic form factor, H., and the 
f 1 

momentum thickness, 8 . The corresponding relation for compressible flow 
was obtained by referencing the parameters in the Ludweig-TilIman skin 
friction relation to the Eckert reference enthalpy. Finally, the normalized 
boundary layer shear distribution. 
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appearing in the moment-of-momentum integral equation, was related to the 
transformed form factor and the incompressible skin friction coefficient 
using the results of Libby, et. al. (Reference 14). The resulting turbu¬ 
lent boundary layer equations, in terms of physical quantities, become 


df 

dx 


= 1.268 


f 

' M 

dM 
_ e 

dx 

1 + (S + 1)H. + A 

x W 1 


e 


J 



(12a) 


and, 


dH. . dM 

_ 1 _ _1_ _e 

dx ‘ " 2M dx 
e 


H. (H. + ir(H. - 1) 


1 + S 


H“: + 4H. - 1 
w (HTTiy^H . +' Y) 


A(Ht - 1) 


H. - 
1 


0.011(H. + 1)(H. - \) c 2 T o 


H. 

i 


'f T 


(12b) 


where, 
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As with the nonlinear laminar boundary layer equation, the quartic 
Runge-Kutta integration scheme is used to solve the system of Equations 
12a and 12b . Unlike the laminar boundary layer method, the heat trans¬ 
fer coefficient is not given by Sasman and Cresci's method. Instead, 
Reynold's analogy is used to relate the Stanton number, St, to 
the skin friction coefficient. 


St =• 


p'u C (T 
e e p r 
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During the current effort, no modifications were made to the form of 

the turbulent boundary layer equations. However, the Mark IV authors again 

replaced T /T , wherever it appeared in the equations, by (p Q /P e ) ' 

used the freestream value of p instead of the local value. The kinematic 

o 


viscosity, 

U 

were corrected. 


was also based on freestream conditions. All such , * iblems 


c• Co mparisons and D i scussion 

Comparisons of the laminar and turbulent integral boundary layer 
methods used in the Mark IV program were made with other simplified method, 
and with experimental data. The comparisons were made to validate the 
FORTRAN coding of the methods, but not to ascertain the validity of the 
methods themselves. However, one obvious deficiency of the equations, as 
codedin the Mark IV program, is that they were derived from the two- 
dimensional boundary layer equations. For arbitrary geometries, it would 
be more appropriate to use the axisymmetric form of the equations which 
accounts, in part, for the spreading of the streamlines. To arrive at the 
axisymmetric form of the equations one must begin with the axisymmetric 
continuity equation. 


a( ouR 

9x 







instead of the two-dimensional form. 


Aku! + = o (i 5 ) 
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used in the derivation of Equations 11 , 12a , and 12b . The variable 

R is the local radius of the axisymmetric body. For a two-dimensional body, 

R ->• ^ and dR/dx = 0, and Equation 14 reduces to the two-dimensional form. 

Utilizing Equation 14 instead of 15 in Sasman and Cresci's analysis, 
it may be shown that Equations 12a and 12b may be expressed in axisym¬ 
metric form by adding the term (- f/R)(dR/dx) to the right-hand side of 
Equation 12a . For arbitrary three-dimensional bodies, however, the inclu¬ 
sion of this term requires the proper choice of the radius R. One possible 
approach is to approximate (i/R) by the surface curvature normal to the 
local streamline direction. The rate of change of this transverse curva¬ 
ture along the streamline is also required. In this manner the problem of 
determining the effect of streamline spreading on skin friction and heat 
transfer is reduced to a geometrical problem. 

In general, the effect of streamline spreading is to increase the 
skin friction and heat transfer coefficients. For example, the laminar 
heat transfer coefficient on a cone is a factor of/3 higher than that on 
a flat plate, for identical inviscid edge conditions and wall temperature. 
Since the laminar and turbulent integral methods used in the Mark IV pro¬ 
gram were derived using the two-dimensional continuity equation, the user 
should be aware that the methods will under-predict the heat transfer on 
axisymmetric or three-dimensional geometries. 

The FORTRAN coding of the integral methods was checked out by first 
using a flat plate geometry. For a zero pressure gradient it may easily 
be shown that Cohen and Reshotko's prediction of the skin friction 
coefficient reduces to 
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The Mark IV skin friction output was in exact agreement with hand calcu¬ 
lations using Equation 16 . Code predictions of heat transfer rate were 
compared with hand calculations using Equation 16 and Reynold's analogy 

q = Pr' 2/3 )p u C (T - T ) (17) 

M y 2 /eep'rw 

where the recovery temperature, T , was based on a recovery factor of 
1/2 r 

r = Pr ' . The code's prediction of q agrees well with the simple Reynold's 
analogy as shown in Figure 10. 

Comparisons of the Mark IV code's turbulent skin friction predictions 
for a flat plate were made with calculations using van Driest's method 
(Reference 15). In Figure 11, the two are shown to be in excellent agree¬ 
ment for T w /T q = 1 (approximately an adiabatic wall), but for nonadiabatic 
walls Sasnidn and Cresci's method seems to predict significantly higher 
values of C f than does van Driest's method. However, the Mark IV code 
output matches calculations presented by Sasman and Cresci in Reference 5. 

The coding of the Mark IV integral methods was also checked out for 
nonzero pressure gradients using a spherical geometry. Laminar heating 
rates at the stagnation point were calculated by the Mark IV program, and 
were compared with calculations using the Fay-Riddel method (Reference 16). 
Although the stagnation point heating rate predicted by the Mark IV code 
actually corresponds to a cylindrical stagnation point, the difference 
between cylindrical and spherical stagnation point heating should only be 
approximately 10%. However, factors of 2 and 3 were observed between the 
Fay-Riddel and the Mark IV predictions, depending on the freestream con¬ 
ditions used. A thorough examination of the Mark IV calculations revealed 
that the error was primarily due to the inability of the Mark IV code to 
accurately predict the stagnation point velocity gradient which strongly 
influences the heating rate. Consistent with the Fay-Riddel calculation, 
the normal velocity gradient computation (local tangent wedge) was re¬ 
placed by the relation 
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Figure 10. Heat Transfer on a Flat Plate - Laminar Boundary Layer 
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where R is the radius at the stagnation point and p^p and p<-p are the 
total pressure and density, respectively behind the normal shock. As 
shown in Figure 12, this approach yields a vast improvement in the stagna¬ 
tion point heat transfer prediction. However, due to funding and time 
limitations this procedure was not made a permanent part of the Mark IV 
program. 
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SECTION IV 


INFORMATION FOR THE PROGRAMMER 

As discussed in Section II the general organization of the modified 
Mark IV program is identical to that of the original version (Reference 1). 
However, some major changes were made to certain components of the program 
including the overlay, which contains the viscous methods, and the overlay 
which traces surface streamlines. Outlined in Subsection 1, Program Struc¬ 
ture, are the organizational changes made to each of the two overlays. 

In Subsection 2, New Local Storage, detailed descriptions of two new 
random access units are given. Unit 50 contains surface spline information 
for each of the geometry Panels and surface property data (P/P^, T/T^, etc.) 
along each of the streamlines. Unit 51 is used to store the coordinates of 
the streamlines for later use by the interactive picture drawing program 
TEKPIC (Reference 2). 

1. Program Structure 

The algorithms developed during the present effort to trace continuous 
surface streamlines required that the subroutines comprising overlay (MARK 
4,2,6) in the original version of the Mark IV program be replaced. Shown 
in Figure 13 is a flow diagram of the new streamline overlay. The names of 
three of the subroutines in the original Mark IV program (STREAM, SFNTRP, 
and VALUE) have been retained in the new streamline overlay since these 
subroutines perform functions which are analogous to those required in the 
original program. However, the actual FORTRAN coding of the new routines 
bears no resemblence to the original coding. Two additional subroutines 
(RNGKTA and TRACE) are included in the new streamline analysis. Given in 
Table 1 are descriptions of all subroutines in the new streamline overlay. 
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Figure 13. Structure of Streamline Overlay 
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Table 1. Overlay (MARK 4,2,6) Subroutine Descriptions 


SUB ROUTINE 

RNGKTA 


SFNTRP 


STREAM 


TRACE 


VALUE 


DESCRIPTION _ 

Subroutine responsible for integrating the streamline 
equations (see Section III - Theory, Equations 5a and 
5b) employs a quartic Runge-Kutta Integration scheme. 
Called from subroutine TRACE. 

Prior to the streamline calculations this subroutine 
maps the two appropriate independent variables [(X,Y), 
(X,Z), or (A,<{))] to the unit u-v plane. The surface 
flow properties (P/P ro , T/T^, M, surface velocity direc¬ 
tion cosines, and the dependent geometry coordinate) 
of each Panel are then fit with the surface spline. 

SFNTRP is called from STREAM one time for each Panel 
of the vehicle. After each call, STREAM places the 
spline information on Unit 50 for later use. 

Main program for the streamline analysis. Responsible 
for initializing the appropriate variables, reading 
user-prepared streamline data (flap, etc.), reading 
surface flow properties from Unit 10, coordinating the 
surface spline calculations, and coordinating the actual 
streamline tracing. Routine STREAM also contains the 
algorithm used to trace "additional" streamlines. 

This subroutine is called one time from the main routine, 
STREAM, for each streamline. Given the starting point 
for a streamline, this subroutine traces it until a true 
origin is reached (DIRECT = -1 on the Streamline Data 
Card, Section II - User's Guide), or until a trailing 
edge is encountered (DIRECT =1). If DIRECT = -1, the 
running lengths along the streamlines are reordered so 
that S = 0 corresponds to the true origin. TRACE also 
retains on Unit 50 the surface properties along the 
streamlines which are subsequently used by the Viscous 
Methods option. 

Given the spline coefficients and the bilinear mapping 
coefficients for a particular Panel, this subroutine 
returns the interpolated values of the surface flow 
properties at a given point. The coordinates of the 
point may be given in the Cartesian system used by the 
Mark IV program, or they may be given in (u,v) coordi¬ 
nates. If the coordinates are specified with respect 
to the Mark IV system, subroutine VALUE also provides 
the corresponding (u,v) coordinates of the point, and 
vice versa. Subroutine VALUE also indicates whether or 
not the given point lies within the boundaries of the 
given Panel. 
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The Viscous Methods option, overlay (MARK 4,2,4), employs two 
aoproaches for viscous related computations: (1) "flat plate" techniques 
applied to a simplified geometry model (Mark III Skin Friction Method), 
and (2) boundary layer calculations made along the inviscid streamlines 
generated by the Streamline Analysis. Outlined in Figure 14 is the inter¬ 
relationship of those subroutines used to compute viscous effects along 
the surface streamlines. Since neither the structure nor the overall 
operation of the Mark III Skin Friction method was modified, those sub¬ 
routines associated with tne Mark III approach are not shown. 

With one exception all subroutines called by INTEG are associated 
only with the integral boundary layer methods which are applied along the 
surface streamlines. Subroutine TEMP, however, employs flat plate skin 
friction methods which may be used along the streamlines in place of the 
integral methods (see ISFM flag on the Boundary Layer Method Control Card, 
Section II - User's Guide). Subroutines CFINPT, SFNTR3, VALU3, and INTI 
u?.ed in the original version of the Mark IV program are not required in 
the modified version and were removed. Otherwise, the subroutine organi¬ 
zation of the Viscous Methods overlay remains unchanged from the original 
structure. 

The two subroutines most affected by the modification of the Viscous 
Methods overlay are INTEG and LAMNAR. Surface flow properties along the 
streamlines are accessed from Unit 50 by subroutine INTEG. (Unit 10 con¬ 
tained the streamline data in the original version of the program.) The 
organization of Unit 50 is discussed in Subsection 2, New Local Storage 
Units. The laminar integral boundary layer equation used previously in 
LAMNAR was replaced with a more general equation, as discussed in Section 
III - Theory. This ordinary differential equation is solved in subroutine 
LAMNAR using the quartic Runge-Kutta method. 

Other changes to the Viscous Methods overlay were minor and do not 
affect the general flow of the program. Such modifications included a 
restructure of the labeled COMMON blocks and the correction of several 
local property calculations (discussed in Section III). 
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Given in Table 2 are brief descriptions of the functions performed by 
the subroutines listed in Figure 14. 


Table 2. Description of Modified Subroutines in Overlay 
(MARK 4,2,4) 


Subroutine 

CURVFT 

FUNCT 

GRADNT 

INTEG 


INT2 

LAMNAR 

LGRNGE 

PRECAL 

PROF IL 


Description _ 

Evaluates a polynomial f(x,y) at a specified point 
(x,y) given the coefficients of the polynomial. 

Contains an expression for the displacement thickness 
6 * in terms of the correlation number n. Subroutine 
FUNCT is used in LAMNAR when the initial correlation 
number must be determined from a user-specified dis¬ 
placement thickness (see also ROOT). 

Computes the gradient of a tabulated function using 
finite difference techniques. 

Subroutine responsible for coordinating boundary layer 
calculations along streamlines. Initializes variables, 
reads user-prepared data and options, accesses the 
surface property data along the streamlines from Unit 50, 
and calls the various boundary layer analyses in a logi¬ 
cal order. 

Yields a dimensionless shape factor based on momentum 
thickness [(0 2 /\> w ) du e /dx]. Subroutine INT2 is used 
in LAMNAR to help locate the transition point. 

Solves the laminar boundary layer equation (see Section 
III -Theory, Equation 11), checks for laminar instability 
and transition to turbulent flow, and computes the initial 
values for the turbulent boundary layer analysis. 

Interpolates a given tabulated function using Lagrange's 
four-point method. 

Performs preliminary calculations for the integral 
boundary layer methods, including the smoothing of data 
along the streamlines, if necessary, and the computa¬ 
tion of all gradients required by the boundary layer 
equations. 

Computes the velocity profiles for the laminar 
(Pohlhausen quartic) and turbulent (power law) boundary 
layers. PROFIL also prints the boundary layer para¬ 
meters computed in subroutines LAMNAR and TURBLN. 
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Table 2 (Cont'd.) Description of Modified Subroutines in 
Overlay (MARK 4,2,4) 


Subroutine 

ROOT 

RUNKUT 

SIMPS 1 

SMTHNA 

SPLINE 

TEMP 

TURBLN 

VISCUS 

For a 


Description ___ 

Locates the root of a given function f(x) within a 
specified interval. 

Solves the coupled differential equations for the 
turbulent boundary layer using a fixed-step quartic 
Runge-Kutta method. 

Integrates a given function f(x) over a specified in¬ 
terval using Simpson's rule. 

Routine used to smooth tabulated data. 

Computes the first and second derivatives of a func¬ 
tion of one variable using a cubic spline technique. 

Primarily responsible for calculating an equilibrium 
wall temperature for the Mark III Skin Friction method. 
Subroutine TEMP may also be used to calculate skin 
friction coefficients along the inviscid surface stream¬ 
lines. 

Coordinates the integration of the turbulent boundary 
layer equations and subsequently calculates other 
boundary parameters of interest (momentum thickness, 
displacement thickness, etc.) 

Main routine for the Viscous Methods overlay. Coordi¬ 
nates the Mark III Skin Friction calculation, but does 
not coordinate viscous calculations along streamlines 
(see INTEG). 


more detailed description of the subroutines see Reference 6. 
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2. New Local Storage Units 

The tracing of continuous surface streamlines over arbitrary vehicle 
configurations requires that a surface interpolation method be available 
for estimating geometrical and surface flow properties at points other than 
the Elements' centroids. As discussed in Section II, Subsection 2 (Surface 
Streamline Tracing), inviscid analyses in the Mark IV program calculate 
surface flow properties only at the Elements' centroids, and the data are 
subsequently placed on random access Unit 10 for use by the streamline 
analysis. The interpolation method of Harder and Desman’s (Reference 3) 
is then used to surface fit the following quantities as a function of two 
appropriate geometrical variables: P/P^, T/T^, M, Vj, V 2 , V 3 , and X3, 
where V-j is the it * 1 component of the unit surface velocity vector and X3 
is the coordinate not used as an independent variable in the surface fit. 

Due to the nature of the interpolation method, each of the seven 
variables are surface fit Panel-by-Panel. Thus, in order to provide a 
complete description of one particular variable's surface distribution 
(e.g., p/p^), a total of N surface fits are required corresponding to the 
N Panels of the vehicle. Therefore, a total of 7 X N surface fits are 
needed to describe the surface distribution of all seven variables over 
the entire vehicle. 

The surface fitting procedure, described in Section III - Subsection 1 
(Streamline Tracing) and in References 1 and 3, is performed in overlay 
(MARK 4,2,6), which contains the streamline analysis, prior to the actual 
tracing of the streamlines. For each Panel of the geometry the coordinates 
of the Elements' centroids and the corresponding flow variables are read 
from Unit 10, fit with the surface spline, and the resulting spline coef¬ 
ficients are saved on random access Unit 50 for later use. The streamline 
calculation begins when the spline coefficients of the last Panel are 
computed and saved. As the streamlines are integrated from one Panel to 
another, a new set of spline coefficients are read from Unit 50 and are 
used to interpolate for the values of the flow variables at any point on 
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the surface of the Panel. Due to the size of the arrays required to hold 
the spline coefficients of each Panel, only the coefficients associated 
with one Panel may be placed in memory at any one time. 

Two types of data are placed on random access Unit 50. Data related 
to the surface splines are stored on the first half on Unit 50, and 
streamline data are placed on the second half of the unit. As shown in 
Table 3, two records are allocated to each Panel for storing the spline 
coefficients and related variables. The Panel information is placed on 
Unit 50 in the same order that the Panel Identification Cards are read 
by program GEOM, the Mark IV geometry program (see Reference 1, Vol. I, 
p. 18). That is, records 1 and 2 correspond to Panel 1, records 3 and 4 
to Panel 2, etc. 


Table 3. Structure of Random Access Unit 50 
(See Table 4 for definitions of all 
variables) 


Record Contents ____ ___ 

1 Ten integer variables related to the surface splines 
of Panel 1. (See first 10 elements of COMMON block 
SURFIT below and Table 4). 

2 Real variables related to the surface splines of 
Panel 1. (See remainder of COMMON block SURFIT below 
and Table 4). 


3 Same integer variables as those in record 1, but as 
applied to Panel 2. 

4 Same real variables as those in record 2, but as 
applied to Panel 2. 

• 

2*NP-1 Same integer variables as those in record 1, but as 

(NP = » Panels) applied to Panel NP. 


2*NP 


Same real variables as those in record 2, but as 
applied to Panel NP. 
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Table 3 (Cont'd.) Structure of Random Access Unit 50 
(See Table 4 for definitions of all 
variables) 


Record 

Contents 


2*NP+1 

ISFLAG(IO) 

Streamline 1 

2*NP+2 

SSTRM(150) 

Streamline 1 

2*NP+3 

XSTRM(150) 

Streamline 1 

2*NP+4 

YSTRM(150) 

Streamline 1 

2*NP+5 

ZSTRM(150) 

Streamline 1 

2*NP+6 

EMSTRM(150) 

Streamline 1 

2*NP+7 

PSTRM(150) 

Streamline 1 

2*NP+8 

TSTRM(150) 

Streamline 1 

2*NP+9 

CFSTRM(150) 

Streamline 1 

2*NP+10 

• 

IPSTRM(150) 

Streamline 1 

• 

t 

2*NP * 1 

ISFLAG(IO) 

Streamline NS 

+10*(NS-1) 

• 



• 

• 

2*NP+1 

IPSTRM(150) 

Streamline NS 


+10*(NS-1) 

+9 

The spline-related data are retrieved from Unit 50 after having first 
placed the FORTRAN names in a single labeled COMMON block: 

COMMON/SURFIT/N1,N2,N3,IU,IW,ID,IORN,ISYM,ICL,IDUM(1), 

1 CPT(4,3),0RGN(3),R0T(3,3),P$IO,THETO,PHIO,UD(500),WD(500), 

2 ULE,WLE,FL0WC(7),FL0WD(500,7),AU,AW,BX,CX,DX,BY,CY,DY,DUM(56) 
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The data for Panel number IPANEL are read from Unit 50 using the 
following statements: 

IREC = 2*1 PANEL -1 

CALL READMS (50, Nl, 10, IREC) 

IREC = IREC + 1 

CALL READMS (50, CPT (1,1), 4600, IREC) 

The spline data for Panel number IPANEL are then available to all sub¬ 
routines in the streamline tracing overlay through the labeled COMMON 
block. 

Following the last record containing spline information are the stream¬ 
line data. Ten records are allocated to each streamline with each record 
containing the values of a particular flow quantity (e.g., T/Tj along that 
streamline. 

Listed below in Table 4 are brief descriptions of the variables stored 
on Unit 50. 


Table 4. Dictionary for Unit 50 Data 
(listed in relative order of 
appearance) 


I. Spline-related Data 

Variable Name Description _ 

Nl Number of centroids in the Panel. 

N2 Number of data points used to generate the surface 

spline (usually, N2 = Nl). 

N3 Number of spline coefficients used to surface fit a 

particular quantity. N3 = N2 + 3, and is the same 
for each of the seven variables that are surface fit. 
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Variable 

IU 

IN 

ID 

IORN 

ISYM 

I CL 

IDUM(l) 

CPT(4,3) 

ORGN(3) 


Table 4 (Cont'd.) Dictionary for Unit 50 Data 
(listed in relative order of 
appearance) 

Name Description _ 

Identifies the first independent variable used in the 
surface fit. Corresponds to the x-coordinate if the 
functional form of the Panel is either Y = F(X,Z) or 
Z = F(X,Y). IU corresponds to the axial component 
A if the Panel's functional form is R = F(A,0), where 
R is a radius and 0 is a circumferential angle. 

IU = 1 always. 

Identifies the second independent variable used in 
the surface fit. IW = 2 if either of the functional 
forms Z = F(X,Y) or R = F(A,0) is used. If the form 
Y = F(X,Z) is used, IW = 3. 

Identifies the coordinate not useo as an independent 
variable in the fit. This third coordinate is fit with 
the surface spline in the same manner that the six flow 
variables are fit. 

Flag indicating whether element data for this Panel 
were input by cross-section or in streamwise strips. 
(See Panel Identification Card description. Reference 
1, Volume I, p. 18). 

Indicates whether or not the vehicle has a plane of 
symmetry. ISYM = 0 if a plane of symmetry exists. 
Otherwise, ISYM = 1. 

Has the value 1 if ISYM = 0 and the vehicle is oriented 
such that the sideslip angle is zero. Otherwise, 

ICL = 0. 

Not used. 

Array containing the coordinates of a Panel's four 
corner points. 

Coordinates of the origin used when the Panel is fit in 
a cylindrical coordinate system, R = F(A,0). If the 
appropriate functional form of the Panel is Y = F(X,Z) 
or Z = F(X,Y) then 0RGN(1), (2), (3) = 0. 
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Table 4 (Cont’d.) Dictionary for Unit 50 Data 
(listed in relative order of 
appearance) 

Variable Name Description _ 

ROT(3,3) Matrix used to transform the Mark IV coordinate system 

to one in which the x-axis is aligned with the axial 
coordinate A in a cylindrical coordinate system. If 
the Panel's functional form is not R = F(A,0) then 
R0T(3,3) is the identity matrix. 

PSIOJHETO Yaw, pitch and roll angles, respectively, used to 

form the matrix R0T(3,3). 

UD(500), The u- and v-coordinates of the centroids of the Panel 

WD(500) (see Section III - Theory). 

ULE,WLE Variables used to describe the orientation of the Panel 

relative to the x-axis of the Mark IV coordinate system. 
Used only when the Panel's functional form is R = F(A,0). 

If the dot product of the unit vectors along the cylindri¬ 
cal axis A and the Mark IV, x-axis is less than .707 (45°), 
then ULE = 0, WLE = -10. Otherwise, ULE = -10., WLE = 0. 

FL0WC(7) Array containing the average value of each of the seven 

variables that are surface fit. 

FL0WD(500,7) Array containing the spline coefficients of each of the 

seven spline-fit variables. 

AU,AW,BX,CX, Certain additive and multiplicative combinations of 

DX,BY,CY,DY the four sides of a Panel. Used in the interpolation 

process. 

DUM(56) Not used. 

II. Streamline Data 

ISFLAG(IO) Array of integer flags. ISFLAG(l) is the number of 

points saved along a given streamline. Remainder of 
array is not used. 

SSTRM(150) Running lengths along a streamline. 

XSTRM(150) The x=, y-, and z-coordinates, respectively, of the 

YSTRM(150) points along a streamline. 

ZSTRM(150) 
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Table 4 (Cont'd.) Dictionary for Unit 50 Data 
(listed in relative order of 
appearance) 


Variable Name 
EMSTRM(150) 

PSTRM(150) 

TSTRM(150) 

CFSTRM(150) 

IPSTRM(150) 


Description _ 

Local Mach number along the streamline. 

P/P^ along the streamline. 

T/T^ along the streamline. 

Not used. 

Panel numbers along the streamlines. 


The organization of Unit 50 does not permit the saving of multiple 
cases. Only the streamline data associated with one set of freestream 
conditions (Mach number, angle of attack, etc.) may be saved at any one 
time. Therefore, if boundary layer calculations are to be made along the 
streamlines, the viscous analysis must immediately follow the streamline 
calculations. 

A second random access device. Unit 51, was added to the Mark IV pro¬ 
gram to provide a communication link with program TEKPIC, an interactive 
computer graphics program designed specifically for displaying Mark IV 
geometries and surface streamlines. If the user-specified variable IST0RE 
is set equal to 1 on the Streamline Data Card (see Section II - Surface 
Streamline Option), then the coordinates of points along the streamlines 
are automatically saved on Unit 51. These streamline data are used only 
by program TEKPIC and do not affect any calculations made in the Mark IV 
program. 

Provisions are made on Unit 51 for storing multiple cases. A maximum 
of 5 streamline distributions may be stored, each corresponding to a 
particular set of freestream conditions, and all bookkeeping is handled 
internally. If the streamline coordinates are to be used subsequently 
by the TEKPIC program, the contents of Unit 51 must be saved on a 
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permanent file device following the execution of the Mark IV program. 
Furthermore, it is not necessary that all 5 streamline cases be placed 
on Unit 51 at once (one execution of the Mark IV program). If, after 
having saved 1, 2, 3, or 4 streamline distributions on one file of a 
permanent file device, it is permissible to attach that permanent file 
device at a later time as Unit 51 and to use the Mark IV program to 
generate new streamline distributions. The new streamline distributions 
are automatically appended to the streamlines already existing on the 
file. However, the streamline data must not be used with program TEKPIC 
until all desired streamline distributions have been placed on Unit 1. 

The structure of Unit 51 is fully documented in Reference 2 for those 
interested in using program TEKPIC. 
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